Quantum catcher - stopping particles of unknown velocities 
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We propose a method to stop particles of unknown velocities by collision with an accelerated wall 
with trajectory ~ y/t. We present classical and quantum mechanical descriptions and numerical 
simulations that show the efficiency of the method. 
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I. INTRODUCTION 



A beam of particles can be slowed down by reflecting 
them from a potential wall (or "mirror") moving in the 
same direction. An early example is the production of an 
ultracoldbeam of neutrons colliding with a moving Ni- 
surface Moving mirrors for cold atom waves have 
been also implemented with a time-modulated, blue- 
detuned evanescent light wave propagating along the sur- 
face of a glass prism [2|, |3[ . More recently, beams of He- 
lium atoms have been slowed down using a Si-crystal on a 
spinning rotor 0,(7]. Also, Rb atoms have been stopped 
with a moving magnetic mirror on a conveyor belt Q, 
which provides a promising mechanism to generate a con- 
tinuous, intense and slow beam of atoms. There is clearly 
a great potential for practical applications of such pro- 
cesses and much interest in their fundamental properties 
and optimization. In most cases the analysis is made 
with classical trajectories, but quantum motion effects 
may become important for ultracold atoms, as shown in 
a recent study on matter- wave/moving- mirror interac- 
tion Q. One more limitation of standard settings so far 
is that the mirror motion is adapted only to stop a spe- 
cific initial beam velocity. We may however consider the 
much more general case in which the initial velocities are 
broadly distributed, or that they are unknown. Stopping 
a beam in these conditions is a much more challenging 
objective, and the central subject of the present paper. 



Can we stop a classical particle of unknown velocity 
with a moving hard wall? The answer is yes and the 
trajectory of the wall is surprisingly simple, as we shall 
see in Sect. |TTJ We shall also explore in Sect. IIIII the 
extent to which the same methodology can be applied in 
the quantum case, and provide numerical examples. 
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II. CLASSICAL ENSEMBLE OF 
NON-INTERACTING PARTICLES 

Let us start by assuming a classical point-particle emit- 
ted at x = 0, t — 0, moving with unknown velocity v > 
along x. A heavy hard wall (compared to the mass of 
the particle) moves also with trajectory x m (t) and initial 
condition x m (Q) =0. If the particle touches the wall at 
time t with velocity v s and the wall is moving with ve- 
locity v m at i, then the final velocity vj of the particle 
after the perfect reflection is 



v f = -v s 



2v r , 



(1) 



This can be easily seen in the reference frame in which 
the mirror is at rest at time t. To fulfill the goal that 
the particle is at rest in the laboratory frame after the 
collision, the mirror velocity should be v m — v s /2 at the 
time of the collision. Since the trajectory of the particle 
with velocity v > is x(t) = vt and the trajectory of the 
mirror is x m (t) , the time t c of the collision is given as a 
solution of x m (t c ) — vt c . The condition for stopping the 
particle is now 



dt 



Using v = x m (t c )/t c we get 



dt 



■He) 



*(*c) 



2t c 



The trajectory for the mirror we are looking for is found 
as the solution of this ordinary differential equation with 
the initial condition x m (0) = 0, 



x m {t) = aVt, 



(2) 



where a > is in principle arbitrary, although its value 
will determine the location and time of the particle- 
mirror collision, which is quite important in a practi- 
cal implementation with limited space and time. Such 
a moving wall will stop particles starting at the origin 
and moving with positive velocity independent of their 
initial velocity, see also Fig. [T] Instead of a we will deal 
with more intuitive parameters: we assume a final time 
t f , this may be a maximal time we are ready to consider 
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FIG. 1: Scheme of the stopping of classical particles: a hard 
wall moving with trajectory x m (t) = d^Jtjtf (solid line); ex- 
amples of two particle trajectories with different initial veloc- 
ities (dashed lines). 




FIG. 2: Scheme of the transformation between initial pa- 
rameters and final parameters, the wall is indicated by the 
black box, the dark gray region is for free motion, and the 
light gray region for motion with collision. The symbols 
provide examples connecting final and initial parameters: 
x s /d = OjVs/vt = — > Xf/d — QjVf/vi = (free motion, 
circles); x s /d — 0,v s /vb = 1 — > xj/d — l,Vf/vb = 1 (free mo- 
tion, boxes); x s /d = 0,v s /v b — 2 — > xj/d = 1/2, Vf/vb = 
(collision, triangles). 



in our experiment. At this time, the wall has moved a 
distance d (see Fig. [T]), so we get a = d/y/tj. Another 
important quantify is Vb := d/tt, the boundary velocity 
for no collision until tf, i.e., a particle with an initial 
velocity v s < Vb will not hit the mirror until tf. 



Let us now examine the effect of the mirror with tra- 
jectory x m (t) = dyjt/tf in more detail, allowing for a 
more general scenario where the initial position of the 
particle is not necessarily zero but a;(0) = x s < 0. We 
are interested in the particle's position xf and velocity 
Vf at tf. It is useful to introduce dimensionless variables 
to simplify the notation in the following, namely 



x v t 

X=-;,v=—,t=—. 
d Vb tf 



By a straightforward calculation we get for the position 
and velocity of the particle at time r = 1 for \s < 0: 



Vf{ X s,Vs) = { 2 Us '" e) 



v s : v 8 > 1 - Xs 
: otherwise 



XfKXs,v s ) - } xs + vs ■ otherwise 



with 



v(Xs, v s ) = — (l + y/l - 4 Xs v s ^ . (4) 

If v s < 1 — Xs the particle keeps its initial velocity, i.e., it moves too slowly to collide with the moving- wall before ti, 
see also Fig. [5] We can also work out the inverse transformation to get the initial position and velocity from the final 
parameters for Xf < 1' 

„ fy. v ,\ = { X (Xf,Vf) ~Vf ■ Xf>v f and v f <0 
sVA -f' i' \ Vf : otherwise 

Xs (Xf,»f) ^ / ^ " ^ + ("/ " :X/>-/and, / <0 



Xf ~ v f 



: otherwise 
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with 



2(x/ ~ f/) 



(5) 



Vf/Vb 




FIG. 3: Final velocity Vf(x s ,v s ), the line at v s /t>{, = 1 —x s /d 
separates the regions with or without a collision before t / . 



The function Vf(x s ,v s ) = VbVf{x s /d 1 v s /vb) = 
VbVf(xs,v s ) is important as it provides the extent to 
which the moving mirror fails to stop the particles when 
they deviate from the ideal conditions \s = and v s > 1. 
This function is shown in Fig. [3j Let us examine the case 
v s > 1 — Xs (and Xs < 0), i.e. a collision has occured be- 
fore t = 1: the final velocity is negative, ^/ < 0; also 



dv s 



= -1- 



< 0, 
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so the absolute value of the final velocity is increasing 
with increasing \x s \ and \v s \. 

Fig. [3] suggests to us a possible strategy to select the 
parameters d and tf and optimize the stopping: assume 
that the initial parameters x s and v s are fixed (they may 
correspond to estimates of the least favorable values ex- 
pected or permitted, such as the farthest distance from 



the origin allowed by the initial geometry of the launch- 
ing conditions and a lower bound for the speed). First 
we choose d and tf such that 



d 

— < Vs 



We have for the final velocity 

Vf(X S ,Vg) 



V h Vf 



X s V s 

d ' v b 



(vb fixed). 



According to Fig. [31 this means that if we make d large 
enough, increasing also tf so that Vb remains constant, 
the final velocity goes to zero, even with imperfect con- 
ditions. The condition Vb(l — x s /d) < v s should be also 
satisfied for the chosen d, so that the particle collides 
with the wall, but this is not a problem because 



Xs 

d 



Vb <C v s (vb fixed). 



In the following we shall discuss the more general case 
in which the initial position and the velocity of the parti- 
cle are characterized by a probability density p s (x s ,v s ). 
The final probability density for position and velocity of 
the particle is given by 



Ps{xf,Vf) =p s [x s {Xf,Xf),V s (v f ,Vf)}. 

In particular we shall consider examples with 

(V-V Q ) 2 (x-X ) 2 



Ps{x,v) 



N 



exp 



2Aw 2 2Air 2 
for x > (JV being a normal- 



for x < 0, and p s (x, v) 
ization constant). 

The initial probability density p s for the parameters 
x a /d = -0.04, Ax/d = 0.008, v /v b = 5.0 and Av/v b 
2.0 is shown in Fig. [4^ (solid line). Also the final prob- 
ability density is shown in Fig. QJi (striped line). The 
integrated probability densities 



Px,s/f(x) = J dvp s/f (x,v) 
Pv,s/f(v) = / dxp s/f {x,v) 



for this example are shown in Fig. UJd and c. Fig. 2J: 
illustrates the slowing down and narrowing of the final 
velocity distribution (dashed line) compared to the initial 
one (solid line). Moreover, a small fraction of particles 
have not hit the wall: they correspond to the final distri- 
bution in the interval < v/vb < 1. One more example is 
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shown in Fig. [5l where a rather drastic stopping can be 




seen. Note that the decrease of the width of the velocity 
distribution results always in an increase of the width of 
the position distribution because the phase-space volume 
is conserved. 
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FIG. 4: Classical setting, x /d = -0.04, Ax/d = 0.008, 
vo/vb = 5.0 and Av/vf, = 2.0; (a) initial probability den- 
sity ps (solid graph) and final probability density pf, (b) p x>a 
(t — 0, scaled by a factor of 1/40) and p x j (t = tj); (c) p VtS 
(t = 0, scaled by a factor of 4) and p x j (t — if). 



III. QUANTUM SETTING 



The dimensional Schrodinger equation for the particle 
and the moving mirror potential is 



d h 2 d 2 

ih—ip(x,t) = - — ^2^( x '*) +V ( x ~ d \/ t / t f ) ^(a;,*)) 



where V(x) is the potential of the wall and the mass 
in the examples below is that of Rubidium (rridim = 
14.19226 • lCT 26 kg). The probability densities do not 
depend on the mass of the particle classically, but in 
the quantum case they do. The initial wavefunction is 
a Gaussian, but not necessarily a minimum-uncertainty- 
product one, 



2(1 + 2iAv 2 fi5) 



[iv (6v - 2x) + 2Av 2 fi{(x - f3) 2 + 2Sv /3)] 



where S = \J AAx 2 - l/(Au 2 ^ 2 )/(2Au), and (3 = x - 
Svq, [i = m/h. The form of Heisenberg's uncertainty 
relation is AxAv > l/(2/j,). 

In the calculations we shall use an infinitely high wall, 
i.e, 

... i f oo : x > 

as well as a more realistic Gaussian wall, with potential 

V G (x) = F exp[-a: 2 /(2A4)]. 

For cold atoms such a potential could be implemented by 
a detuned laser. 



The first example corresponds to the parameters 

x a = —2 /im, 

Ax = 0.4 /im, 

v = 3.125cm/s, 

Av = 1.25 cm/s. 

We choose a boundary velocity Vb = 0.625cm/s vq and 
d = 50/^m, so that tf = 8 ms. We have xo/d = 0.04 <C 1. 
The height of the Gaussian potential is set to Vo/Ti = 
3.75 • 10 6 /s and its width is Axy = 0.4/im, both parame- 
ters must be chosen such that the particle is reflected with 
high probability. The initial and final quantum probabil- 
ity densities are shown in Fig. Qualitatively, we see 
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FIG. 5: Classical setting, x /d = -0.003, Ax/d = 0.0003, 
Vo/vb = 3.0 and Av/vb = 1.24; (a) p x , s (t — 0, scaled by a 
factor of 1/500) and p x j (t = tf); (b) p v , a (t = 0, scaled by a 
factor of 50) and p x j (t = tf). 
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FIG. 6: Quantum setting, see text for parameters; (a) p x , s 
(t — 0, dashed line, scaled by a factor of 1/40), p x j with an 
ideal wall (t — tf, thick solid line), p x .f with a Gaussian wall 
(t — tf, thick dotted line); (b) p V)S (t = 0: scaled by a factor 
of 4, dashed line), p v j with an ideal wall (t = tf, thick solid 
line), p v ,f with a Gaussian wall (t = tf, thick dotted line) and 
P*,f (t = !)• 



the same result in the quantum case as in the classical 
case, with quantum interference fringes superimposed. 

Another spectacular example can be seen in Fig. 
which shows the initial and final quantum probability 
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FIG. 7: Quantum setting, see text for parameters; (a) p XlS 
(t — 0, dashed line, scaled by a factor of 1/100), p x j with an 
ideal wall (t = tf, thick solid line), p x j with a Gaussian wall 
(t = tf, thick dotted line); (b) p v<s (t = 0: scaled by a factor 
of 100, dashed line), p v j with an ideal wall (t = tf, thick 
solid line), p v j with a Gaussian wall (t — tf, thick dotted 
line) and p x j (t — 1). 

densities for 



X = 


."! fjxa, 


Ax = 


0.3 /xm, 


vo = 


3 mm/s, 


Av = 


1.24 mm/s 


d = 


1000/im, 




Is, 


v b 


1 mm/s. 



We have v b < v and x /d = 0.003 < 1. The height of 
the Gaussian potential is set to VqJU = 3 ■ 10 4 /s and its 
width to Axv — 0.8 //m. 

IV. SUMMARY 

We have proposed a method to stop particles of un- 
known velocities or ensembles of particles with initial 
velocity spread by reflecting them from an accelerated 
wall with trajectory ~ \ft. For classical particles, the 
stopping is perfect if they are emitted from the origin at 
t = 0. If the particle starting point is not at the ori- 
gin or if it is too slow, it is not perfectly stopped in a 
finite time. Explicit expressions are given for the final 
velocity; they also show how to mitigate, and even sup- 
press in a limit, the effect of non- ideal initial conditions. 
We may expect a similar behavior for a quantum wave 
packet and indeed, using numerical simulations with re- 
alistic and experimentally accessible parameters, we have 
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illustrated the efficiency of the method and discussed its 
bounds. 
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